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New measurement of the ^Li/^Li isotopic ratio in the extra-solar 
planet host star HD82943 and line blending in the Li 6708 A 

region. 
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Abstract. The presence of possible blends in the spectral region of the Li resonance line at 6708 A in solar-type 
metal-rich stars is investigated using high resolution and high signal-to-noise spectroscopic observations. Our 
analysis does not confirm the identification of a weak absorption feature at 6708.025 A with the low excitation 
Ti I line proposed by Reddy et al. (2002) . Our spectrum synthesis suggests that the unidentified absorption is 
most probably produced by a high excitation Si I line originally proposed by MuUer et al. (1975). Reanalysis of the 
''Li/^Li isotopic ratio in HD 82943 was performed by taking the Si I line into account and using new VLT/UVES 
spectra of HD 82943 with a signal-to-noise ratio close to 1000. We confirm the presence of ^Li in the star's 
atmosphere while the updated value for the isotopic ratio is /(^Li) = 0.05 ± 0.02. 
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1. Introduction 

An interesting opportunity for testing the planet (s) inges- 
tion scenario is offered by a ^Li test (Israelian et al. l2U01|l . 
This approach is based on looking for an element that 
should not appear in the atmosphere of a normal star, 
but would be present in a star that has swallowed a planet. 
We have proposed recently that a light isotope, ^Li, is an 
excellent tracer of any planetary matter accretion. The 
detection of ^Li in HD 82943 by our group (Israelian et 
al. I2001|l was considered as convincing observational evi- 
dence that stars may accrete planetary material, or even 
entire planets, during their main sequence (MS) lifetime. 
Other explanations (such as stellar flares or surface sports) 
of this phenomenon were ruled out (Israelian et al. l2UUl|l . 
It has been suggested (Sandquist et al. .2002j) that ^Li can 
be used to distinguish between different giant planet for- 
mation theories. 

Nuclear reactions destroy the lithium isotopes (^Li and 
^Li) in stellar interiors at temperatures 2 x 10^ (^Li) and 
2.5 X 10^ K (''Li). Convection cleans the upper atmosphere 
of Li nuclei by transporting them to deeper, hotter layers 
where they are rapidly destroyed. Young low mass stars 
are entirely convective and most of the primordial Li nuclei 
are burned in their interiors in a mere few million years. 
However, many solar-type stars preserve a large fraction of 
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their initial atmospheric ^Li nuclei. According to standard 
models (Forestini 1994. Montalban & Rebolo 2002J , at 
a given metallicity there is a mass range, where ^Li but 
not ^Li is being destroyed. These models predict that no 
^Li can survive pre-MS mixing in metal-rich solar-type 
stars. Observations seem to support this scenario since no 
clear detections of ^Li have been reported in solar-type 
solar-metallicity stars (Andersen, Gustaffson & Lambert 
119841 Rebolo et al. ll986|l . We note, however, that the mass 
and the depth of the convection zone also depend on the 
metal content of the star, and for this reason several old 
metal-poor stars have preserved a fraction of their initial 
^Li nuclei (Cayrel et al. lT^ Hobbs, Thorburn & RebuU 
[l999) . This suggests that any detected ^Li in a metal- rich 
solar-type star would most probably signal an external 
source for this fragile isotope. 

Observations of ^Li are extremely difficult for several 
reasons. First of all, it is a weak component of a blend of 
the much stronger doublet of ^Li with an isotopic sepa- 
ration 0.16 A. In the case of metal-poor halo stars with 
[Fc/H]< —1, blending of the Li line with other weak ab- 
sorptions is not expected, and also the placement of con- 
tinuum does not pose serious problems. ^Li has been un- 
successfully sought in many metal-poor stars but unam- 
biguously detected in only few (Hobbs et al. 119991 Cayrel 
et al.IEll Nissen et al. ITO^ . The methods used in the 
analysis of ^Li have been widely discussed in the litera- 
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ture (e.g. Nissen et al. [TM^ Cayrel et al. [TUM Hobbs 
et al. 11999(1 . In metal-rich stars the identification of any 
possible weak blends in the region of the Li absorption 
becomes crucial (King et al. I1997f) . Recently Reddy et 
al. (2002') updated the line list of Lambert, Smith & 
Heath (1989) based on the analysis of the solar spectrum 
of Kurucz et al. l|TMHl and Hinkle et al. iPTOjl . They 
claimed that a previously noticed weak absorption in the 
solar spectrum at 6708.025 A (Miiller, Peytremann & De 
La R,eza ll97^ belongs to Til. With this assumption their 
analysis of the Li feature in HD 82943 did not confirm the 
presence of ^Li. 

With the goal of establishing the nature of the absorp- 
tion at 6708.025 A and better understanding the line list 
near Li, we have obtained high resolution, high S/N spec- 
tra for several metal-rich stars with no detectable Li line. 
Our targets cover a wide range of effective temperature, 
allowing a more reliable identification of the lines in the 
spectral region of the Li feature. We have also obtained 
new spectra of HD 82943. Here we present new analysis 
of the ^Li/'^Li ratio in HD 82943 and critically discuss the 
line list of the Li region. Our observations do not con- 
firm the Ti I line of Rcddy et al. (20(21) and support our 
previous claim for detection of ^Li. 

2. Observations 

We used the 8.2 m VLT Kueyen (ESO, La Silla), the 
UVES spectrograph and an EEV detector to obtain spec- 
tra of HD 82943 during the night of 2001 March 5. The 
0.3" slit yielded a resolving power of i? 105 000 as mea- 
sured from the FWHM of Th-Ar lamp lines. 

The 4.2 m William Herschel Telescope on La Palma 
(Canary Islands, Spain) and the UES spectrograph were 
used to observe several solar-type stars with large Li 
depletions in their atmospheres. Our first targets were 
HD 171888 (F8 V), HD4747 (G8 V), HD 217580 (K4 V), 
HD 217107 (G3 V) and HD 210277 (G4 V). These stars 
were observed on 2000 July 8. The CCD SITEl (2148 x 
2148 pixel), grating E31, slit width 1" (providing a resolv- 
ing power of 55 000) and the central wavelength of 5140 A 
were used during this run. Our spectra had an S/N ratio in 
the range 300-350, which was sufficient to guarantee the 
detection of any spectral lines with a minimum equivalent 
width of ~1.5 mA. These data were used by us (Israelian 
et al. 2001) in order to check the presence of blends near 
6708 A with the EW larger than 1.5 mA. More targets 
(Table 1) have been observed in order to reduce this limit 
to below 1 mA. The same configuration was used in 2001 
October to observe HD 16141, HD 22049, HD 75732A and 
HD 217107. In the next observing run (2002 June 20), we 
used a 4096 x 2048-pixel EEV CCD with a pixel size of 
15 /im, the E31 grating, a slit width of 1" and a central 
wavelength of 7215 A. Ten exposures of HD 82943 were 
taken on 2001 February 5 with a mosaic of two EEVs, the 
E79 grating and a slit width of 0.7" in order to obtain a 
high quality of spectrum of this star. 



The SARG high resolution echelle spectrograph at the 
3.5 m Telescopio Nazionale Galileo (La Palma) was used 
during two runs in 2001 August and 2002 May. The spec- 
tra were obtained with the yellow grism and span the 
wavelength range 4600-7820 A at a resolving power of 
-57000. The CCD was a mosaic of two 4096 x 2048-pixel 
EEVs with a pixel size of 15 /im. 

AU our VLT/UVES, WHT/UES and TNG/SARG 
spectra were reduced using standard IRAF^ routines. 
Normalized flats created for each observing night were 
used to correct the pixel-to-pixel variations and a Th-Ar 
lamp was used to find a dispersion solution. 

Some observations were made with the FEROS spec- 
trograph at the ESO 1.52 m telescope in La Silla. The 
mosaic of two 4096 x 2048-pixel EEV CCDs were used 
to observe HD 1461 and HD13 445 on 2000 November 10 
and HD 192263 on 2001 October 31. The spectra were flat- 
fielded, calibrated with a Th-Ar lamp and reduced using 
MIDAS routines. 

3. Spectrum synthesis and stellar parameters 

In this article we used models of atmospheres provided by 
Kurucz (1992) and the spectrum synthesis code MOOG 
(Sneden 1973). We have compared abundances of various 
elements computed with our atmospheric models and with 
those used by Reddy et al. (1999) (i.e. Kurucz 1995). We 
found that differences induced in the abundances by the 
two types of models are less than 0.01 dex and can there- 
fore be neglected. Solar abundances of chemical elements 
were taken from Anders & Grevesse (|1989|l . We have also 
made an extensive use of the VALD database (Kupka 
et al. 1999). The parameters of some of the planet host 
stars listed in Table 1 were slightly updated with respect 
to the values given in Santos et al. (2001, 2002). Three 
stars in our sample (HD 1461, HD 198084 and HD 200790) 
are not known to harbour planets and were observed be- 
cause of the absence of Li. We have attempted to es- 
timate the parameters of a cool subgiant, HD 137759, 
from the set of Fe lines used in our previous analysis 
(Santos et al. r2DDT)) . However, given the low temperature 
of this object, our final values had large error bars and we 
have decided to use the parameters based on colours and 
accurate parallax measurement from Hipparcos (Mallik 
1999). Two Li-poor dwarfs, HD 198084 and HD 200790, 
were taken from Boesgaard et al. (|2001fl . The parame- 
ters of HD 200790 were estimated in the same way as for 
the other stars in our sample. HD 198084 is a metal-rich 
F8IV-V (Boesgaard et al. 2001, Malagnini et al. I^nnn|l 
spectroscopic binary (Griffin .1990jl . The iron lines in our 
spectrum were strongly blended owing to the presence of 
the second star and we did not make a new estimate of 
the stellar parameters. 

^ IRAF is distributed by National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation, USA. 
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Table 1. Observing log and stellar parameters adopted from the references listed in the last column: (1) Santos 
et al. (2002), (2) Santos et al. ^UU^, (3) Santos et al. (2003), (4) MaUik ifTM^ . (5) this article, (6) Boesgaard et 
al. l|MIT|l . The S/N is calculated near 6706.5 A. 
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4. The line list in the Li region 

Accurate measurements of the wavelengths and oscilla- 
tor strengths of the Li lines are available in the literature 
(Sansonetti et al. 1995). Unidentified blends in the Li re- 
gion of the Sun and other cool stars have been discussed by 
several authors (Miiller, Pcytrcmann & De La Reza'1975 
Lambert, Smith & Heath 1989, Andersen et al. .l984.; King 
et al. nWzjl . Miiller et al. ^WJ^ were the first to notice 
the weak feature at 6708.025 A in the solar spectrum and 
attributed it to the Sii line with loggf = —2.7, x = 6 eV. 
The line list of Lambert et al. 1)1989(1 has been extensively 
discussed and modified by King et al. 1)1997(1 . Following 
Miiller et al. (1975), these authors used a fictitious high 
excitation Si i line with xi =6 eV in order to account for 
the weak absorption in the spectrum synthesis. 

The line was also considered by Nissen et al. (1999) 
in their analysis of metal-poor disc stars. According to 
Nissen et al. (1999), the equivalent width (EW) of the 
6708.025 A absorption is 0.6 mA in the solar spectrum. 
Reddy et al.'s 1)2002)1 identification of this absorption with 
the Ti Hue (loggf = -2.252, x = 1-88) imphes a shghtly 
different value EW = 0.75 mA. While this feature is too 
weak to affect the determination of the Li isotopic ratio 
in metal-poor halo stars (where it certainly disappears), 
it may cause problems in solar metallicity or more metal- 
rich stars. Apart from this unidentified feature, there are 
other weak lines of Cri 6707.64 A, Cei 6707.74 A and 
CN 6707.816 A with uncertain gf values. It is interesting 
that the best synthetic spectra fits to the Li region of 
different stars can be achieved only by altering the gf 
values of these weak lines by different amounts. This led 



King et al. 1)1997)1 to propose that various uncertainties in 
the analysis of the solar spectrum allow Li isotopic ratios 
as high as 0.1. 

Most of the lines in the list of Reddy et al. 1)2002)1 come 
from the original compilation of Lambert et al. 1)1989)1 . 
However, the Til 6707.752 A line was not listed in any 
of the references cited by Reddy et al. 1)2002)1 . On the 
other hand, the Vi 6707.563 A line in the list of Reddy 
et al. iPTU^ appears under 6707.518 A in the VALD- 
2 database (Kupka et al. 1999). The gf values of many 
lines have been modified by Reddy et al. 1)2002)1 in order 
to fit the solar spectra of Kurucz et al. (1984), Hinkle et 
al. (2000) and 16 Cyg B. In most cases there is a clear 
disagreement with the list of King et al. l)1997)l . who em- 
ployed the same spectrum synthesis code, model atmo- 
spheres and, as a matter of fact, the same targets (i.e. the 
Sun and 16 Cyg B). Despite this. King et al. 1)1997)1 found 
for CN 6707.816 A loggf = -1.74 while Reddy et al. 
(2 002) arrived at loggf = —2.317. Futhermore, King et 
al. p997)l have considerably increased the loggf value of 
a Ceil line at 6707.74 A and obtained log 5/ = -1.86 
while Reddy et al. 1)2002)1 decreased it to —0.8 dex af- 
ter Lambert et al. (1993) and obtained a value loggf = 
-3.81. The Vi 6708.07 A line with loggf = -2.33 from 
Lambert et al. l|T^89,) appears under Vi 6708.094 A with 
loggf = -3.113 in the list of Reddy et al. ^U^. The 
latter took it from the list available in the web site of R. 
Kurucz. The gf values of the Cri 6707.64 A line in King 
et al. (1997) and Reddy et al. (2002) differ by 0.14 dex. 
The source of these discrepances was not investigated by 
Reddy et al. (i2002i . 
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Our own synthesis of the solar spectrum suggests gf 
values close to those proposed by King et al. ^997). 
For example, when doing a synthesis of ^Li in HD 82943 
(Israelian et al. 2001) we set log gf = -1.81 for the CN 
6707.816 A line, which is similar to loggf = —1.74, as 
found by King et al. (1997). In addition to these weak 
lines, there is a Til line at 6707.964 A with x = 1-88 
eV and loggf = —6.903, which appears in the VALD-2 
database. This fact probably led Reddy et al. H2002|l to 
assume that the feature at 6708.025 A also belongs to Ti 
and has a similar excitation energy but much larger os- 
cillator strength. For some reason Reddy et al. (2002) did 
not incorporate the Ti line at 6707.964 A in their final list. 
The two Til lines at 6708.025 and 6708.125 A of Reddy et 
al. H2002|l have never been tested in stars other than the 
Sun and 16 Cyg B. 

It is clear that the Li line list in the solar spectrum 
is not well calibrated. Different authors obtain different 
adjusted gf values even when they use the same tools and 
spectra. Our goal is not to make a better fit to the solar Li 
region since this does not guarantee a unique identification 
of the spectral lines. Instead, we have decided to use the 
same list as Reddy et al. H2002|l and see if we can explain 
observations of stars other than the Sun and 16 Cyg B. 

5. The Li region in the Sun and in Arcturus 

The unidentified line at 6708.025 A is clearly seen in the 
solar spectrum (Fig. 1). In the same figure we show a 
small spectral window from McDonald's atlas of a K III 
giant Arcturus observed with a resolving power 150000 
and S/N - 1000 (Hinkle et aL lMHljl . Most of the absorp- 
tions in this small region appear in the spectra of both 
stars, which have very different temperature, gravity and 
metallicity. It is interesting that the feature at 6708.025 
A does not become stronger in Arcturus. It may still ex- 
ist in the spectrum but is severely blended with a strong 
absorption centred at 6708.094 A. 

The line at 6708.094 A is certainly not a redshifted 
Til 6708.025 A of Reddy et al. ij^UH^jl . It also cannot be 
one of the TiO lines observed in much cooler stars (Luck 
I1977|l and taking into account in our synthesis of three cool 
giants. The wavelength of this line corresponds to the V i 
at 6708.094 A in the VALD-2 database and also appears 
in the list of Reddy et al. 1)2002(1 but with an oscillator 
strength that does not fit observations of Arcturus and 
other cool giants. Both wings of this line are disturbed 
by blends. The figure also demonstrates that blending is 
not the only problem in these studies. The exact location 
of the stellar continuum becomes a critical issue when we 
deal with absorption features as small as 1 mA. 

We know that echelle spectra are not calibrated in the 
absolute flux scale, and that, therefore, the placement of 
local continuum levels is either visually estimated after 
fitting the extracted spectra with polynomials, or results 
from a detailed spectral synthesis. We have marked sev- 
eral points (C1-C6) in Fig. 1. which indicate very narrow 
windows theoretically free from spectral lines. These win- 
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Fig. 2. Observed (filled small squares) and synthetic spec- 
tra of three metal-rich stars: HD 217107, HD 75732A and 
HD 145675. The Ti abundances were computed for [Ti/H] 
= (thin continuous), 0.2 (dashed), 0.4 (thick continu- 
ous) in HD 217107, [Ti/H] = (thin continuous), 0.15 
(dashed), 0.3 (thick continuous) in HD 75732A and [Ti/H] 
~ (thin continuous), 0.2 (dashed), 0.4 (thick continuous) 
in HD 145675. 
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Fig. 3. Observed (filled small squares) and synthetic spec- 
tra of the Arcturus, HD 177830 and HD 137759. The 
Ti abundances were computed for [Ti/H] = —0.2 (thin 
continuous), —0.3 (dashed), —0.5 (thick continuous) in 
Arcturus, [Ti/H] — —0.2 (thin continuous), 0.0 (dashed), 
0.3 (thick continuous) in HD 177830 and [Ti/H] = -0.4 
(thin continuous), 0.0 (dashed), 0.4 (thick continuous) in 
HD 137759. 

dows are used in metal-poor stars to make a fine tuning 
of the continuum and estimate the S/N. Points CI, C2, 
C5 and C6 mark the solar continuum fitted by Kurucz et 
al. (1984). We do not know how the continuum was placed 
in Arcturus but it appears from the figure that only C5 
and C6 can be considered as continuum points. Neither is 
it clear whether C2 and C4 are continuum points in the so- 
lar spectrum. The point C3 lies on the continuum in those 
stars that do not have CN and Li absorptions at 6707.78 
A. The strength of the unknown line at 6708.025 A and 
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Fig. 1. High resolution and high S/N spectra of the Sun and Arcturus in the Li region. The unidentified feature is 
centred on 6708.025 A in the Sun but appears at 6708.094 A in the spectrum of Arcturus. 

Table 2. Abundance of Ti (column 4) as derived from the EWs of two Tin lines (columns 2 and 3) measured in our 
spectra. The EW of the unidentified feature (column 5) at 6708.025 A was computed assuming that it is due to the Ti 
line of Reddy et al. ()20()2p with abundances from column 4. Observed EWs of unidentified line are listed in column 5 
with 3 a errors estimated using a Cayrel formulae fCavrei ri988|l . Values in parentheses in column 6 indicate that the 
line is not centred on 6708.025±0.01 A. The last two columns provide the Si abundance in our stars from Bodaghee et 
al. (2003) and the EW of the unidentified feature assuming that it belongs to Sii with x — ^ and log gf = —2.97. 
It is not clear if the 6708.025 A feature in HD 200790 is real. 
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the Vi absorption at 6708.28 A define whether or not C4 
is placed on the local continuum. The region between C2 
and the Si i line at 6707.01 A is depressed in the solar spec- 
trum as well as in the spectra of many other stars probably 
owing to some weak, still unknown, blends that start to 
appear in cool and very metal-rich stars (e.g. HD 145675 in 
Fig. 2). Depending on the exact location of the points Cl- 
C6 one may derive different EW values for the unidentified 
feature at 6708.025 A. However, as we will show below, 
even this uncertainty does not smear the huge differences 
between observations and spectrum synthesis when the Ti 
hnc of Reddy et al. H2002|l is considered. 

6. The spectral feature at 6708.025 A 

6.1. The Ti I line of Reddy et al. (2002) 

Let us investigate how the Til 6708.025 A line of Reddy 
et al. H2002|l fits our observations. We have computed Ti 
abundance in the sample of metal- rich stars (Table 1) us- 
ing two lines of the ionized Ti from Santos et al. (2000). 
The EWs of these lines and derived abundances are listed 
in Table 2. The difference in abundance from the two lines 
was always less than 0.03 dex. Our tests with the neutral 
Til lines from Santos et al. (2000) showed that the abun- 
dances derived from Til and Tin agree within 0.1 dex 
in F8-G2 stars while there is a systematic difference at 
lower temperatures. The abundances obtained from the 
Tin lines are about 0.1-0.15 dex smaller than those de- 
rived using Til. Because of a small non-LTE effect on 
Til, the use of Tin hues is preferred (N. Shchukina 2002, 
personal communication). In addition, the abundances de- 
rived from Ti n lines are insensitive to the errors in Tcff 
(Santos et al. 2000). In any case, even if we suppose that 
our Ti abundances are uncertain by up to 0.3 dex (which 
is very unlikely, given the errors considered in this work), 
this will not change our final conclusion regarding the re- 
liability of the line proposed by Reddy et al. H2002|l . 

Using the abundances from the Ti n lines, we have es- 
timated the EW of the feature at 6708.025 A, assuming 
that it belongs to the line proposed by Reddy et al. I|2()02(l . 
Predicted and observed EWs of the Ti line are listed in 
Table 2. We see an obvious disagreement at low Tcff: the 
predicted line is much stronger than the observed one. In 
Fig. 2 we show spectrum synthesis for three metal-rich 
stars with different parameters. The abundances of V, C, 
Si and other elements were scaled with Fe and then mod- 
ified by small amounts (always less than 0.1 dex) in order 
to fit the observations. HD 217107, with Toff — 5655 K, is 
more similar to the Sun than the two other stars and ap- 
parently for this reason the agreement between predicted 
and observed profiles of 6708.025 A is good. There is small 
room for making adjustments in continuum level in order 
to accommodate [Ti/H] — 0.3 and fit the observations. 
This plot alone cannot be used to argue against the iden- 
tification proposed by Reddy et al. 1)2002(1 . As a matter of 
fact, the observed and predicted EWs of the Ti line agree 
within errors in HD 217107 (Table 2). The discrepancy is 
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Fig. 4. The unidentified feature at 6708.025 A has EW < 
1 mA in HD 198084 and HD 200790. 



larger in cooler objects HD 145675 and HD 75732A. We 
also note that the centre of the Ti line is slightly shifted 
to the blue in the spectrum of HD 145675. This shift and 
the disagreement between predicted and observed EW be- 
come much larger when we study K-type stars (Fig. 3). 
One of our targets is the well-known red giant Arcturus. 
The atmospheric parameters and the abundances of ele- 
ments in this star were taken from the literature (Peterson 
et al. 1993). This figure demonstrates that the predicted 
Ti line becomes strong and its wavelength clearly does 
not correspond to the absorption observed at 6708.094 
A. Although our computations rule out the [Ti/H] ratios 
—0.2 and —0.4 (thin continuous lines in Fig. 3) considered 
by us as the lower limits in the synthesis of HD 177830 
and HD 137759, we provided them in order to show the 
wavelength shift between the fictitious Ti line and the ab- 
sorption at 6708.094 A. These plots clearly demonstrate 
that the feature at 6708.025 A disappears at low temper- 
atures while a different line at 6708.094 A (not clearly 
observed in the solar spectrum) becomes stronger. The 
wavelength of this line corresponds to V i in the VALD-2 
database. Two lines at 6708.025 and 6708.094 A appear as 
a single absorption feature in stars with Tcs = 5000-5400 
K. It should be possible to separate them in very high 
resolution and high S/N spectra of stars with intrinsically 
narrow absorption lines. 

When estimating the EW of the Ti line in HD 82943 
we found another disagreement with Reddy et al. (2002)1. 
who state that the combined strength of the two Ti i lines 
(6708.025 and 6708.125 A) in the spectrum of HD 82943 
is 2.2 mA, while our analysis using exactly the same tools 
and parameters shows 1.1 mA, in which the 6708.125 A 
line contributes only 0.2 mA. The same disagreement ap- 
pears when we consider the solar spectrum. Using the list 
of Reddy et al. ^UU^, we found that the EWs of the Til 
lines at 6708.025 and 6708.125 A are 0.75 and 0.17 mA, re- 
spectively. Thus, the combined EW — 0.92 mA is at odds 
with the 1.5 mA claimed by Reddy et al. 1)2002)1 . 
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Fig. 5. The spectrum of HD 82943 (small filled squares) 
and two synthetic spectra computed without Li. The first 
synthesis (solid line) is computed without CN, and as- 
suming [Si/H] = 0.1 and [V/H] — 0.3. Computations with 
[C/H] = [N/H] = 0.35, [Si/H] = 0.27 and [V/H] = 0.8 are 
presented with a dashed line. The position of the ^Li line 
is indicated. 

6.2. Identifying the line 

We have attempted a new identification of the 6708.025 
A feature knowing that the strength of the line does not 
increase towards cooler temperatures. It is clear that any 
molecular or a low excitation line from a neutral metal will 
behave like the Ti line of Reddy et al. ( 2002 ) . Since this is 
not observed, we are left with two possibilities. Either the 
transition belongs to an ion, or it is produced by a neutral 
element with a high excitation energy of the lower level. 

The most abundant ions at these temperatures belong 
to Fe, Si, Mg, Ca and Ti. Reliable experimental line lists 
exist for all of these ions and are available in VALD-2. 
Most of the Si ii linesin the Li region have very high exci- 
tation energies of between 10-16 eV. The same is true for 
Mgii (x = 11-13 eV) and Can (x = 7-9.7 eV). Tin has 
several transitions with x = 1-9-3.1 eV. A Tin line with 
X — 1.9 eV and log gf — —4.1 will reproduce the solar fea- 
ture with EW = 0.75 mA. This line will also fit the spectra 
of the coolest objects in our sample, such as HD 137759 
and HD 177830, where it will be weaker than 2 mA. In 
HD 82943 it wiU have EW = 1.1 mA, which is not really 
different from the strength of the neutral Ti line (0.9 mA). 
However, small discrepancies may appear in hotter stars. 
For example, given the Ti abundance in HD 200790, we 
anticipate about 1 mA absorption which we think is not 
observed in the high S/N spectrum (Fig. 4). The marginal 
feature observed at 6708.028 A is about 0.6 mA and it is 
hard to judge by eye if it is real. On the upper panel 
of the same graph we show the spectrum of HD 198084. 
Unfortunately, this star is in a binary system and we can- 
not make any strong statements regarding the strength of 
the unidentified absorption. However, we would still ex- 
pect to see the blend of two unidentified lines at 6708.025 
A. In fact, there is a broad absorption around 6708.0 A, 
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most probably caused by the shifted unidentified lines but 
its EW does not exceed 1 mA. 

We have studied a spectrum of the hot Li-poor star 
HR 7697 (Teff=6588 K, logg = 4,24, [Fe/H] = 0.01) with 
S/N ^ 550 kindly provided by Ann Boesgaard (Boesgaard 
et al. 2001). The ionized line of Ti in this object would 
produce about 1 mA absorption, which we could not find 
in the data. These tests make the idea of an ionized Ti 
line less attractive, although, given the errors in the data 
and our analysis, we still cannot completely rule out this 
possibility. We also note that it will be much harder to 
identify the line if it belongs to an ion of an clement whose 
abundance does not scale with Fe. 

Rare Earth elements were also consid- 
ered. The best candidates are perhaps two 
Cen lines Hsted in the DREAM database 
( http://www.umh. ac.be~astro/dream.shtml) with 
wavelengths 6708.077 and 6708.099 A. The fines are 
weak (log gf = —2.57 and log gf — —2.12) and have 
excitation energies xi = 2.25 and 0.7 eV, respectively. 
Let us suppose that the 6708.077 A fine is responsible for 
the feature at 6708.025 A. The 6708.077 A absorption 
is absolutely negligible in the solar spectrum, and in 
order to produce 0.7 mA it must have log gf — 0.43. 
This very large correction is not allowed by accurate 
radiative lifetime measurements (Zhang et al. 2001). If we 
still suppose that the adjustment is possible, these two 
lines will produce about 4 mA absorption in Arcturus. 
However, the synthetic spectra will not reproduce the 
observed profile at 6708.094 A unless we modify the 
wavelengths. In brief, we can hardly force these two lines 
of Ce n to explain observations. 

Our synthetic spectra also demonstrate that the gf 
values of Vi fines at 6708.094 and 6708.280 A in the 
list of Reddy et al. ((2002)1 do not explain observations 
of HD 137759, HD 177830 and Arcturus (Fig. 3) assum- 
ing that the abundance of V scales with Fe (Sadakane 
et al. 2002). Our tests show that even adopting [V/Fe] 
= —0.5 in these stars, we still cannot fit the feature at 
6708.280 A. Such low abundances of V are excluded from 
the analysis of other spectral lines. On the other hand, the 
fit to the same line in HD 82943 was achieved by boosting 
the abundance to [V/H] = 0.8. This high value is ruled 
out from a synthesis of different Vi lines. We faced simi- 
lar problem for the line at 6708.094 A but in the opposite 
direction. The newly added Til line at 6708.125 A in the 
list of Reddy et al. H2002|l cannot salvage the situation 
as it appears too weak in our targets. We also note that 
the V I line at 6708.280 A does not appear in the VALD-2 
database (Kupka et al. 1999). 

6.3. The Si I line of M Li Her et al. (1975) 

Given the large number of high excitation Si lines in the Li 
region (see VALD-2 database), we have decided to test if 
the Si fine at 6708.025 A proposed by Miiller et al. ifTHTHIl 
can explain our observations. In fact, the high excitation 
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Table 3. Mean abundance [Si/H] = 0.27 ± 0.02 was computed from 11 Si lines observed in S/N-1000 VLT/UVES 
and WHT/UES spectra of HD 82943. 
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5665.56 


4.92 


— 1.980 
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0.27 
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4.93 
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—2.020 
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0.30 
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0.28 
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57.3 


0.26 


5948.55 


5.08 


-1.110 


109.7 


0.26 


6125.02 


5.61 


-1.520 


46.1 


0.26 


6142.49 


5.62 


-1.480 


50.9 


0.29 


6145.02 


5.61 


-1.400 


54.9 


0.25 


6155.15 


5.62 


-0.750 


113.0 


0.27 


6721.86 


5.86 


-1.090 


65.4 


0.26 



Si line that appears at 6707.01 A (Fig. 1) is very similar 
to the one proposed by Miiller et al. (|1975|l . It behaves 
exactly in the same way in stars with different tempera- 
tures as the absorption at 6708.025 A. Given the EW {— 
0.75 mA) of the 6708.025 A feature in the solar spectrum, 
we have estimated log gf — —2.97 for the Sii line with xi 
= 6 eV proposed by Miiller et al. Ijl975|l . Our gf value is 
slightly different from that of Miiller et al. 1)1975(1 owing to 
the differences in the model atmospheres and the adopted 
abundance of Si (Anders & Grevesse ll989|l . which is 0.05 
larger than the value used by Miiller et al. (1975) . We 
have computed the Si abundance in our stars using the 
lines from Bodaghee et al. (2003) listed in Table 3. The 
new list is an extended and improved version of the one 
considered by Santos et al. H2000|l . The abundance errors 
were always less than 0.07 dex and do not affect our anal- 
ysis. Predicted EWs of the Si line are listed in the last 
column. As we can see, there is good agreement between 
observations and synthesis. 

In conclusion, we think it is more likely that the 
spectral feature at 6708.025 A belongs to a high exci- 
tation transition of some neutral atom. The small spec- 
tral window near 6708 A contains many high excitation Si 
lines; therefore, the identification proposed by Miiller et 
al. 1)1975(1 seem to be the best choice. However, we cannot 
rule out the possibility that the 6708.025 A line belongs to 
an ion or to a rare element. More high quality observations 
are necessary to tackle this problem. 

7. The ^Li/^Li ratio in HD 82943 

HD 82943 was newly observed at VLT/UVES with a re- 
solving power 105 000 and a signal-to-noise ratio close 
to 1000. Reanalysis of the ^Li/^Li ratio has been car- 
ried out using the same tools as in Israelian et al. (2001) 
and the line list of Reddy et al. 1)2002(1 . The only change 
we have made in this list is a replacement of the Ti i 
line at 6708.025 A by the high excitation Sii line with 
\og gf = —2.97 and x/ = 6 eV. This change does not help 
in giving a higher ^Li/ ''Li ratio in HD 82943 since the EW 
of the Si I line in HD 82943 is even larger than that of the 



Ti of Reddy et al.lMF{see Table 2). Furthermore, an ac- 
curate abundance of Si in HD 82943 needs to be calculated. 
For this purpose we have used eleven Si lines from Table 

3 and computed their EWs using very high S/N spectra 
obtained with VLT/UVES and WHT/UES. In Table 3 we 
list the lines, measurements and resulting abundances. We 
found [Si/H] = 0.27 ± 0.02 in excellent agreement with the 
value given by Sadakane et al. (2002) after accounting for 
differences in stellar parameters. 

The identification of the chemical clement is also im- 
portant, although Reddy et al. (2002j have ignored this 
point by stating that the identification of the pair as Ti i 
is not critical to the analysis. The EWs of the Ti line 
of Reddy et al. 1(^00^ are 0.45 and 0.9 mA in HD 82943 
for [Ti/H] = and [Ti/H] = 0.3, respectively. This fac- 
tor of 2 difference in the EWs is equal to 0.05 correction 
in ^Li/''Li ratio if the total EW of the Li line is 10 mA. 
This is a significant difference, and given the fact that in 
many stars [Ti/H] can be larger than 0.3, the abundance 
correction may affect the ^Li/^Li ratios by at least 0.03 
dex, depending on the strength of the ^Li line. 

To demonstrate the blending of the Si and CN lines in 
the Li feature of HD 82943, we have computed synthetic 
spectra without considering any of the Li isotopes. Figure 
5 shows that CN and Si make a significant contribution to 
the total EW of the Li line. Given the high C and N con- 
tent of HD 82943 (Santos et aL OUUni Sadakane et al. 00021 
one may consider the abundances of these elements as free 
parameters in the analysis. The strength of the Si line 
may affect the synthesis as well. It is also important to 
realize that this line cannot exactly mimic the ®Li compo- 
nent at 6708.0728 A because their wavelengths differ by 
0.05 A. 

In obtaining the broadening parameters we used the 
same Fe i lines as in our first paper (Israelian et al. 12001(1 . 
From the analysis of two Fe lines (Fig. 6) we ob- 
tained a mean wsini = 1.6 ± 0.05 km s^^ and Vmac — 

4 ± 0.3kms~^, in good agreement with our previous re- 
sults and also with those of Reddy et al. 1(2002(1 . The syn- 
thetic spectra were also convolved with a Gaussian func- 
tion representing the instrumental profile. Given the low 
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Fig. 6. Observed and computed profiles of two Fe lines 
and results of analysis. Filled dots in the upper panel 
are observations and continuous lines are synthetic spectra 
derived from the analysis. Dashed and solid lines in the 
lower panel are the analysis of Fei 6703.6 and 6705.1 
A lines, respectively. 

rotational velocity in HD 82943, it is possible to obtain fits 
using simple Gaussian functions. As we have noticed be- 
fore (Israelian et al. 120011 see also Reddy et al. 12002(1 . 
the final values of ®Li/^Li are not affected by using a 
pure Gaussian broadening function. Average broadening 
parameters obtained from the Fe lines were used to derive 
the ^UfU ratio. 

The free parameters in a standard analysis of the 
^Li/^Li ratio in HD 82943 were the total abundance of 
Li, the wavelength shift, the ^Li fraction and the position 
of the continuum. The final ®Li/''Li ratio is not strongly 
affected by the blue CN absorption, and after few a it- 
erations we fixed the CN abundance equal to [C/H] = 
[N/H] — 0.35. Small adjustments of the continuum level 
not exceeding 0.2% (given the S/N of the spectrum) were 
allowed to optimize the fits. The abundance of the Si line 
was also considered as a free parameter in the range [Si/H] 
= 0.27 ± 0.1 but after a few iterations it was fixed to 
[Si/H] = 0.27. Additional computations were carried out 
allowing for large variations of v sin i and Vmac exceed- 
ing the 3(7 limit derived from the Fe lines. We found that 
these effects do not change our final results. The minimum 
X^ was derived by optimizing all free parameters for each 



value. The new analysis using Reddy et al.'s H2002|l line 
hst and the Si hnc yields /(^Li) = 0.05 ± 0.02 (Fig. 8, left 
panel). The best fit for /(^Li) — has an rms of 2.869 x 
10~^, which is significantly larger than the rms (9.91 x 
10-4) of the best fit with /(^Li) = 0.05. Wavelength off- 
sets in the best fits for /(^Li) = and /(^Li) = 0.05 were 



Fig. 7. Comparison of the observed (filled large dots) and 
synthetic spectra of HD 82943 corresponding to ^Li/''Li = 
0.05 (continuous) and ®Li/^Li ~ (dotted) isotopic ratios. 
They correspond to the best fits with /(^Li) = and a 
wavelength offset of — 0.65 km s"^ (small dots) and to the 
f(^Li) = 0.05 with a wavelength offset of — 0.44kms~^ 
(continuous line). Fits to the blue wing of the Li profile 
can be improved if we adopt the wavelengths of CN lines 
from Brault & Miiller (1975). The CN hues observed in 
the arc spectrum by these authors appear at 6707.55 A 
while Reddy et al. (2002) and Lambert et al. (1993) list 
them between 6707.464 and 6707.529 A. The residuals (O- 
C) of the observations after subtraction of the synthetic 
spectra are shown. 
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Fig. 8. Results from the analysis for the ^Li/^Li ra- 
tio. Ax^ = l,4,and 9 correspond to 1,2 and 3a confidence 
limit is (dotted lines). Continuum adjustments were al- 
lowed within 0.2 % (left panel) and 0.5% (right panel). 



— 0.65 km s ^and— 0.44kms ^ , respectively. These values 
agree perfectly with those reported by Cayrel et al. (1999) 
in their ^Uf'U analysis of the halo star HD 84937. The 
synthetic and the observed profiles for two cases are shown 
in Fig. 7. Given the strength of the Li hne (about 50 mA), 




these shifts are consistent with observations of Fe Hnes in 
the solar spectrum (AUcnde Prieto & Garcia Lopez ni)98|l . 
Since the scatter of the velocity shifts of Fe lines in the 
solar spectrum is large (ranging from —0.8 to 0.2kms~^), 
one cannot make any predictions for the blueshifts of the 
Li lines based on a statistical trend observed for the so- 
lar Fe lines. So the best option still remains to keep the 
wavelength shift of the Li line as a free parameter. 

The minimum is not affected by large variations 
in the continuum either. Relaxing the upper limit of con- 
tinuum change from 0.2 to 0.5 % does not change our 
conclusion (Fig. 8, right panel). The final value of /(^Li) 
is also not sensitive to the Si abundance. If we repeat the 
analysis fixing the continuum but changing the abundance 
of Si, we will find that the best fit to the observed profile 
is achieved for [Si/H] = 0.27 (Fig. 9, right panel). The 
same analysis for /(^Li) = yields [Si/H] = 0.51 (Fig. 9, 
left panel). Such a high abundance of Si is ruled out by 
our data. However, even if we assume that it is possible, 
the quahty of the fit when /(^Li) = 0.05 and [Si/H] = 
0.27 is about 4 times better (i.e. is smaller) compared 
with the case when /(^Li) = 0.0 and [Si/H] = 0.51. On 
the other hand, the Si i line may suffer from convective 
blueshifts as well. To account for this effect, we have re- 
peated the analysis considering the wavelength of the 
Si line as another free parameter. We found that the effect 
is very small. For example, when /(^Li) = 0.05, the —8 
mA wavelength shift will introduce a negligible change in 
X^ of about 5%. This would also imply a small increase in 
the Si abundance equal to [Si/H] = 0.31. These changes 
do not affect our final results. 

Reddy et al. I)2()()2|l noted in their paper that the Ti 
lines contribute 2.2 mA (while our own estimate is EW = 
1.1 mA) into the 50 mA Li line, which, according to these 
authors, corresponds to ^Li/^Li = 0.06. Given that with- 
out this blend the ratio was ^Li/^Li = 0.12 (Israelian et 
al. 120011 . one should wonder what is causing the remain- 
ing 0.06 contribution. If the contribution of the blend is 
only 0.06, then we are still left with another ^Li/^Li — 
0.06. The blended line is not strong enough to completely 
eliminate the contribution from ^Li. Although Reddy et 
al. H2()02(l have stated that our spectra agree very well, 
they did not explore this puzzle. 
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Fig. 10. Upper panel: comparison of the observed (filled 
squares) and synthetic spectra of HD 82943 and a Gen A. 
Lower panel: comparison of the observed (filled squares, 
Kurucz et al. 1984) and synthetic spectra of the Sun using 
the list of Reddy et al. with the wavelengths of GN lines 
from Brault & Miiller (1975) and the Si line discussed 
in this article. There are missing lines at 6707.2 with a 
total EW < 0.3 mA (see also Reddy et al. 2002; Nissen et 
al. 1999; Miiller et al. 1975), which are clearly seen for the 
first time in the spectrum of a very metal rich cool star 
HD 145675 (Fig. 2). These weak lines do not affect our 
analysis (i.e. the continuum placement) as they disappear 
in stars with Toff > 5900 K. 



In order to judge the validity of our spectral synthe- 
sis and the placement of the stellar continuum, we have 
synthesized spectra of a Gen A and HD 82943 in the 7 A 
region near the Li doublet (Fig. 10). Observations of a Gen 
A with a resolving power of 105.000 are taken from Rebolo 
et al. (1986) and have S/N >500 in the Li region. As stellar 
parameters for a Gen A we have adopted Tcff — 5800 K, 
logg = 4.31, = 1.1 kms-\ [Fe/H] = 0.24 (Ghmielewski 
et al. 1992) and abundances were taken from the litera- 
ture (Abia et al. 1988; Santos et al. 2000; Ghmielewski et 
al. 1992; King et al. 1997). 
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Given that we are not able to make a unique identifi- 
cation of the chemical element responsible for the blend, 
the ^Li/^Li found from the analysis cannot not be con- 
sidered to be a final value. The abundance of Si usually 
scales with Fe in planet host stars (see for example Santos 
et at. I^nnni Gonzalez et al. EUIIT! Sadakane et al. 2002, 
Takeda et al. 12001(1 . However, if the spectral line belongs 
to an element which is under- or overabundant with re- 
spect to Fe, then the ^Li/^Li ratio will be different. In most 
cases the elemental abundances are not very different from 
that of Fe (typically within ±0.1 dex) and therefore the 
overall effect on the ^Li/^Li ratio will not be large. 

8. Conclusions 

High quality observations of metal-rich stars with unde- 
tectable Li have essential consequences, first for the iden- 
tification of weak lines in the Li 6708 A region, and sec- 
ond for the analysis of the ^Li/^Li isotopic ratio. Our re- 
sults clearly rule out the Til line at 6708.025 A suggested 
by Reddy et al. I)2002|l . Alternatively, we propose that 
this feature can be attributed to the Si i line of Miiller et 
al. (1975). The identification and the presence of another 
Til hue at 6708.125 A in Reddy et al. (2002) line list is 
also not justified and cannot be verified with our dataset 
as this line is extremely weak in all our targets. We also 
find that two hues of Vi at 6708.094 and 6708.280 A from 
the same list have incorrect parameters. However, these 
Vi lines and the Til line at 6708.125 A cannot strongly 
affect the ^Li/^Li ratio in stars with T^s > 5500 K. 

Our analysis of the new VLT/UVES spectra 
(S/N=1000), taking into consideration the Sii line at 
6708.025, gives an isotopic ratio ^UfU = 0.05±0.02 in 
HD 82943. This implies a contribution to the total ab- 
sorption in the Li region of 1.3 mA and 2.7 mA due to Si i 
and ^Li, respectively. If no Si i blend is considered in the 
analysis we obtain ^Li/^Li=0.08, which is ~ 33 % smaller 
than the value obtained in our first analysis (Israelian et 
al. 2001) based on a spectrum with S/N=500. Since ^Li 
can be construed as a fossil record of rocky matter ac- 
cretion during planetary system formation, further high 
quality observations of this fragile isotope may offer us 
a unique and invaluable insight into the past dynamical 
history of extra-solar planetary systems. 
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